Relativistic positioning: four-dimensional numerical 
approach in Minkowski space-time. 
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Abstract We simulate the satellite constellations of 
two Global Navigation Satellite Systems: Galileo (EU) 
and GPS (USA). Satellite motions are described in 
the Schwarzschild space-time produced by an ideal- 
ized spherically symmetric non rotating Earth. The 
trajectories are then circumferences centered at the 
same point as Earth. Photon motions are described in 
Minkowski space-time, where there is a well known rela- 
tion, Coll, Ferrando & Morales-Lladosa (2010), between 
the emission and inertial coordinates of any event. 
Here, this relation is implemented in a numerical code, 
which is tested and applied. The first application is 
a detailed numerical four-dimensional analysis of the 
so-called emission coordinate region and co-region. In 
a second application, a GPS (Galileo) satellite is con- 
sidered as the receiver and its emission coordinates are 
given by four Galileo (GPS) satellites. The bifurcation 
problem (double localization) in the positioning of the 
receiver satellite is then pointed out and discussed in 
detail. 

Keywords relativistic positioning systems; methods: 
numerical; reference systems 



1 Introduction 

Nowadays, in order to design an operative relativistic 
positioning system (RPS), we need: (i) inertial coor- 
dinates (x^ ) labeling events in an appropri- 
ate reference system, (ii) four satellites, whose posi- 
tions are known at any time, which broadcast their 
proper times by means of electromagnetic signals, (iii) 
detectors which receive the proper times from the four 
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satellites at the same time. These times are the emis- 
sion coordinates of the reception event, and (iv) the 
transformation from emission to inertial coordinates, 
which localizes the reception event in the inertial sys- 
tem of reference (relativistic positioning). In any RPS, 
some 4-tuples of emission coordinates may be received 
in two different positions at two distinct coordinate 
times; namely, each of these 4-tuples leads to two real 
and different reception events in Minkowski space-time. 
If one of these 4-tuples of proper times is received by 
an observer, new information -apart from the emission 
coordinates- is required to choose one of the two pos- 
sible locations. The study of this problem, in realistic 
cases, is one of the main goals of this paper 

Although current positioning systems are based on 
Newtonian physics, relativistic post-Newtonian correc- 
tions are performed if necessary; however, RPS should 
be based on relativistic principles from the beginning. 
In the proposed schemes, the proper times of four satel- 
lites (emission coordinates) are sent by means of elec- 
tromagnetic signals to the receiver, whose inertial coor- 
dinates can be found -from the emission ones- by using 
fully relativistic equations. Realistic four-dimensional 
(4D) implementations of the transformation derived by 
Coll, Ferrando & Morales-Lladosa (2010) require nu- 
merical calculations. See Puchades & Saez (2011) as 
a preliminary numerical application of this transforma- 
tion. Here, numerical codes based on the same trans- 
formation have been designed and tested. Results ob- 
tained with them are described in next sections. The 
transformation from emission to inertial coordinates 
(Coll, Ferrando & Morales-Lladosa 2010) in Minkowski 
space-time -where it is asumed that the electromag- 
netic signals propagate- uses the position of the four 
satellites when they emitted their proper times. The 
circular motion of GPS and Galileo satellites has been 
simulated by using the Schwarzschild space-time cre- 
ated by an idealized spherically symmetric Earth (see 
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next section); thus, the satellite positions may be cal- 
culated at any given time with a good enough accu- 
racy. Of course, the formalism must be extended to gen- 
eral relativity (Cadez & Kostic 2005; Bini et al. 2008; 
Teyssandier & Le Poncin-Lafittc 2008; Cadcz, Kostic & 
Delva 2010; Bunandar, Cavcny & Matzner 2011; Dclva, 
Kostic & Cadez 2011) to include gravitational fields, ac- 
celerated frames, and so on, in both satellite and photon 
motions. Although a fully rclativistic scheme should 
work without any reference to an Earth based coordi- 
nate system (Tarantola et al. 2009), realizations of this 
type of scheme have not been implemented up to now. 

In the study presented here, Earth and other possible 
obstacles to light propagation are not considered at all 
and, moreover, signals broadcast by the satellites may 
be detected at any distance by ideal receivers which 
have not a threshold for detection. In this way, only 
the space-time structure is taken into account in our 
discussion. Finally, a spherically symmetric Earth is as- 
sumed and the satellites describe circular orbits whose 
centers must coincide with that of the Earth. These 
orbits are covered as it must be done in Schwarzschild 
space-time. Afterward, our study may be generalized 
to include obstacles, realistic detectors and non sym- 
metric distributions of mass inside Earth. 

A coherent terminology for relativistic positioning 
may be found in Coll, Ferrando & Morales-Lladosa 
(2010), where concepts as emission coordinates, grid, 
emission coordinate region and co-region, and so on, 
were rigorously defined. Discussions in some previ- 
ous papers (CoU 2001; Bahdcr 2001; Coll 2003; Rovclh 
2002; Bahder 2003; Coll & Pozo 2006a; Coh, Ferrando 
& Morales 2006b, c; Coll, Ferrando & Morales-Lladosa 
2010, 2011) led to synthesize the mentioned terminol- 
ogy and to define the foundations of relativistic posi- 
tioning. A sketch of the emission region in 3D (sec- 
tions of constant coordinate time) is presented in Coll, 
Ferrando & Morales-Lladosa (2012). It is based on a 
system of three static satellites symmetrically placed in 
the vertexes of an equilateral triangle. These satellites 
broadcast their proper times (emission coordinates) to 
receivers which live in a Minkowskian space-time with 
one time dimension and two space ones. Realistic 4D 
studies of emission coordinate regions and co-regions, 
with four dynamical satellites, have not been performed 
yet. The first of these studies is developed here. 

Quantities G, Mq, t, and t stand for the gravitation 
constant, the Earth mass, the coordinate time, and the 
proper time, respectively. Greek (Latin) indexes run 
from to 3 (1 to 3). Quantities rjafj are the covariant 
components of the Minkowski metric tensor. Units are 
defined in such a way that the speed of light is c = 1. 

The paper is organized as follows, in Sect. 2, the mo- 
tion of the GPS and Galileo Global Navigation Satel- 



lite Systems (GNSSs) is simulated, the method used 
to obtain the emission (inertial) coordinates from the 
inertial (emission) ones in Minkowski space-time is de- 
scribed in Sect. 3 (4). In the first (second) of these 
sections we use analytical (numerical) techniques. Sec- 
tion 5 contains numerical results. The emission coor- 
dinate regions (co-regions) defined in Coll, Ferrando & 
Morales-Lladosa (2010) are described in Sect. 5.1 (5.2) 
for various sets of four satellites, and the possibility of 
finding the position of a satellite by using other four 
satellites -as emitters- is discussed in Sect. 5.3. Fi- 
nally, Sect. 6 contains a general discussion and some 
comments about perspectives. 

2 The satellites 

Two GNSSs are considered: GPS in USA and Galileo 

(under construction) in the EU. They provide the spa- 
tial coordinates and the universal time of any event on 
Earth. Position coordinates are calculated thanks to 
information received from satellites into orbit around 
Earth. The GPS constellation has = 24 satellites 
and it is arranged in six different orbital planes (four 
satellites per plane), each of them inclined an angle 
ctin = 55 degrees with respect to the equator. To ob- 
tain around two orbits per day, the satellites are placed 
at an altitude h = 20200 Km. The Galileo constel- 
lation is composed by 27 satellites (n^ = 27), located 
in three equally spaced orbital planes (9 equally spaced 
satellites in each plane). The inclination of these planes 
is ain = 56 degrees and the altitude of the circular or- 
bits is /i = 23222 Km; thus, the orbital period is about 
14h. 

Satellites move in Schwarzschild space-time to take 
into account Earth gravity. The trajectory of any satel- 
lite is assumed to be a circumference of radius R, which 
has the same center as Earth. In Schwarzschild space- 
time, the angular velocity on these circumferences is 

= (GM^/i?^)^/^, so in a coordinate system attached 
to Earth center, the coordinates of a given satellite A 
may be written as follows: 

x\ = R [cos a A (t) cos + sin a a (t) sin V' cos 9] 
Xa = — i? [cos ttyi (r) sin i/> — sin q:^(t) cos V' cos ^] 
.T^ = —R sin a A (t) sin 

x\ = 7r . (1) 

The factor 7 calculated up to first order in GM^/R is 
given by the relation 



3 



and the angle 

aA(r) = a AO - ^Jt (3) 

localizes the satellite on its trajectory. Finally, 9 and 
tp are Euler angles associated to two systems of spa- 
tial axis: the axis {x^,x'^,x^) trivially associated to the 
standard angular Schwarzschild coordinates, and a sec- 
ond set of axis {x'^,x''^,x'^), which is chosen in such a 
way that (a;'^, cc'^) coincides with the orbital plane con- 
taining the trajectory of the satellite under considera- 
tion. Angle 6 = 2?: — Qj„ is the same for all the satellites 
of a given GNSS, whereas angle V takes on the values 
tjj = {2ir/nso){j — 1), where Uso is the number of satel- 
lites per orbital plane and the natural number j labels 
these planes. Evidently, angle tp is the same for all the 
satellites of a given orbital plane. For any satellite, an- 
gles 9, tjj and aAo arc constant. The last angle defines 
the position of satellite A a.t t = =0. This angle 
may be arbitrarily chosen for a satellite in each orbital 
plane and, then, the remaining aAo angles may be fixed 
in such a way that all the satellites are equally spaced 
on their common trajectory. 

3 BVom emission to inertial coordinates 

Events of interest are always simultaneous observations 
of four satellites. The inertial coordinates of one of 
these events are denoted x = {x^ , x'^ , x^ , x'^) = {x,t). 
The emission coordinates arc the four proper times, t^, 
codified in the satellite signals, where index A numer- 
ates the satellites. 

The coordinates {x^,x'^,x^) of the satellite A, at 
emission time t^, are denoted ja- Since the world lines 
of the satellites are known, quantities 7a may be calcu- 
lated for arbitrary proper times. Vectors = 7a — 74 
(with index a running from 1 to 3) define the rela- 
tive position between satellite a and satellite 4, which 
is hereafter the emitter of reference. The numeration 
of the satellites and, consequently, the choice of the 
fourth satellite are arbitrary. We may say that vectors 
Ca define the internal satellite configuration at emis- 
sion times. There arc incrtial coordinates associated to 
times T^, if and only if, the so-called emission- reception 
conditions, Coll, Ferrando & Morales-Lladosa (2010), 
are satisfied. These conditions may be written as fol- 
lows: 

r?„;3e«ef>0, r,„^(e« - e?)(ef - e^) > , (4) 

for any value of indexes a and b which run from 1 to 
3. If these conditions are satisfied, we may look for the 
inertial coordinates. 



In Minkowski space-time, the general transformation 

from emission to incrtial coordinates was derived in 
Coll, Ferrando & Morales-Lladosa (2010); it may be 
written as follows: 

2 

a; = 74 f*^ = , (5) 

(y* • x) + Wiy* ■ x)^ - yix'^ 

where vectors x and y» may be calculated from ei, 62, 
and 63 (internal satellite configuration). The configu- 
ration vector X = *{ei A 62 A 63) (dual of a double exte- 
rior product) is orthogonal to the hyperplane containing 
the four ja emission events. Vector = (^, i?)/(^ ■ x), 
where H) stands for the interior product, may be cal- 
culated from any arbitrary vector ^ satisfying the condi- 
tion^'X ^ and from the bivector iJ = [{eA-eA) /'2\E^ , 
where E'^ = *(e2 A 63), E'^ = *(e3 A ei), and E^ = 
*(ei A 62). 

The above transformation is the solution of a system 
of four equations (hereafter the main system). Each 
equation expresses that the distance from ja to x van- 
ishes; so two types of solutions appear. The first type 
corresponds to signals emitted from the satellites at 
times and simultaneously received at position x and 
time t (emission or past-like solutions). The second 
type describes a signal emitted from position x and 
time t and received by the satellites at times (re- 
ception or future-like solutions). In any RPS we are 
only interested in the first type. Hereafter, solutions of 
this type are also called positioning solutions. 

For x^ 7^ 0, there are two sets of inertial coordi- 
nates corresponding to e = -|-1 and e = — 1. Moreover, 
for < 0, only one of the two sets of inertial coor- 
dinates corresponds to a positioning solution. In the 
case > 0, the number of positioning solutions may 
be either two or zero, in the first case, there are two 
different receptors (located at different places), which 
would receive the same four emission times from the 
same satellites. In the second case, there are two future- 
like solutions (zero positioning solutions). Finally, for 
X^ = there is only a positioning solution correspond- 
ing to e = -|-1. 

In case x^ < 0, the positioning solution satisfies the 
condition — ^ < for any A, whereas the inequali- 
ties tA —t > are satisfied for the future-like solution. 
Since the satellite world lines are known, the inertial 
coordinate tA may be calculated at any proper time 
and, consequently, the sign of tA — t may be used to 
identify the positioning solution. 

If one has a bifurcation problem with x^ > and two 
positioning solutions xi and X2 (Schmidt 1972; Abel 
& Chaffee 1991; Chaffee & Abel 1994; Grafarend & 
Shan 1996), the receiver should have a criterion to se- 
lect its true inertial coordinates. One of these criteria 
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-proposed in Coll, Ferrando & Morales-Lladosa (2011); 
Coll. Ferrando & Morales-Lladosa (2012)- is as follows: 
consider a conical surface with the receiver at the ver- 
tex which contains three of the four satellites and then, 
take either one sign of e or the opposite one depending 
upon whether the fourth satellite is inside or outside 
the cone, respectively: of course, the receiver should 
have devices to measure angles. Since the coordi- 
nate times of the two positioning solutions are different, 
the receiver may select the true positioning solution by 
using a clock. Only the coordinate time of the true solu- 
tion will be identical (close enough taking into account 
the clock accuracy and possible positioning errors) to 
that given by the receiver clock. 

4 From inertial to emission coordinates 

Given the inertial coordinates, x" , of an event in 
Minkowski space-time, its emission coordinates, r^, 
may be numerically calculated. Let us now describe 
the method we have implemented to perform this cal- 
culation. Since emission and reception events must be 
on a null geodesic, we can write the following algebraic 
equations 

V<.pK-xI{t^)][xP-x^^{t^)]=Q, (6) 

in which the proper times are the unknowns. The 
solution of these equations are the emission coordinates 
T^. This solution may be easily obtained by using 
the Newton-Raphson method (Press et al. 1999) plus 
Eqs. (1) and (3). 

After obtaining the four emission coordinates r"^, we 
can use Eq. (5) to recover the inertial coordinates we 
had initially chosen. We use multiple precision in the 
code designed to solve Eq. (6) and also in the numerical 
calculations based on Eq. (5). If a precision of forty dig- 
its is required, we have verified that the parameter fix- 
ing the precision of the Newton-Raphson code may be 
adjusted to recover 39 digits after computing the initial 
inertial coordinates with Eq. (5). This test ensures that 
our numerical calculations leads to very accurate emis- 
sion (inertial) coordinates starting from inertial (emis- 
sion) ones. In other words, we have very accurate codes 
to calculate inertial coordinates from (5), as well as to 
get emission coordinates by solving Eq. (6). The sec- 
ond of these calculations -based on Newton-Raphson 
method- is more time consuming. 

5 Numerical results 

In this section, the emission region and the co-region 
-defined in Coll, Ferrando & Morales-Lladosa (2010)- 



are numerically studied for the first time, in the case of 
realistic satellite configm-ations (GPS and Galileo). 

Numerical methods and codes based on the transfor- 
mations between inertial and emission coordinates have 
been described in Sees. 3 and 4. 

The space 7?.^ containing all the 4-tuples of proper 
times is called grid. Given an arbitrary point of the 
grid {t^ ,t'^ ,t'^) , the question is: would these proper 
times be received in some point of Minkowski space- 
time? In other words, are there positioning solutions of 
the main system for these proper times? If affirmative, 
the points in Minkowski space-time (receivers) would 
belong to the so-called emission region and the chosen 
grid point to the emission co-region. 

The emission region and its co-region are 4D spaces 
and, consequently, graphic descriptions require the 
study of a set of appropriate 3D sections. In the 
emission region (co-region) the chosen 3D sections are 
characterized by the condition a;^ = constant (r^ = 
constant). 

A few words about the figures are worthwhile to en- 
sure a right intuitive interpretation. 

In this paper, the HEALPIx [hierarchical equal area 
isolatitude pixelisation of the sphere) package (Gorski, 
Hivon & Wandelt 1999) is used to depict the figures. 
This pixelisation was designed to construct and analyze 
maps of the cosmic microwave background. It is useful 
to display any quantity depending on the observation 
direction (pixel). The number of pixels is 12 x N^i^eJ 
where the free parameter A^side takes on even natural 
values. In Figs. 1 and 2 corresponding to the emis- 
sion region (5.1), we will take Ngide = 16 (3072 pixels), 
whereas in the Figs. 3 to 5 of the co-region (5.2), we 
will use Nside = 32 (12288 pixels). The angular area 
of any pixel is ~ 13.43 (~ 3.36) squared degrees for 
Nside = 16 {Nside = 32). In the case of Nside = 16 
inside = 32) the pixel is close to sixty four (sixteen) 
times the mean angular area of the full moon, but its 
shape is not always the same (as it is appreciate in the 
figures). Pixels are more elongated in the polar zones. 

Finally, the pixelised sphere is shown by using the 
moUwide projection, in which, the frontal hemisphere 
is represented in the central part of the figure. The 
opposite semi-sphere is divided in two zones which are 
displayed in the lateral parts of the same figure, whose 
edges represent the same back semi-meridian. 

5.1 Emission region structure 

The emission region of four satellites is the zone of the 

space-time where proper times from them (emission co- 
ordinates) may be received. In Minkowski space-time, 
any point, whose emission coordinates (see section 4) 
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satisfy the relation > for the four satelhtes, be- 
longs to the emission region; evidently, is the time 
at which the satellite A started to emit. Signals from 
the four satellites will reach any position {x^,x^,x^) - 
whatever its distance to Earth may be- for a certain 
value of the coordinate time. The larger the distance 
to Earth, the greater this time. 

Two types of reception events are distributed in the 
emission region. The first type is characterized by the 
condition ^ and, consequently, there is only a po- 
sitioning solution (hereafter, single positioning); how- 
ever, the second type corresponds to > and, in 
such a case, there are pairs of events corresponding to 
the same emission coordinates (hereafter, double posi- 
tioning or bifurcated location). An additional criterion 
is necessary for positioning (see section 3). We are inter- 
ested in the distribution of these types of events inside 
the emission region. 

In order to perform an exhaustive study of various 
emission regions we proceed as follows: (1) a reception 
event is selected. It occurs at coordinate time tn in a 
point on the Earth surface with coordinates (x^, a;^, x^). 
From the coordinates {xl,Xg,xl,tii), the emission co- 
ordinates and the quantity ni^y be calculated (see 
Sees. 3 and 4). Since the relation < is always 
satisfied on Earth, the selected event corresponds to a 
single positioning. In all the cases considered along the 
paper, coordinates {xl always the same, and 

they correspond to a point on Earth with the spheri- 
cal coordinates 6e = 60° and = 30°. Three 4-tuples 
of satellites (labeled 1, 2, and 3) are studied. Finally, 
various times tji covering a period of the Galileo satel- 
lites are considered for 4-tuple 1; (2) the hypersurface 
x'^ = tn (emission region section) is studied. In order 
to do that, the HEALPIx pixelisation is used to define 
a set of 3072 directions (see Sect. 5). A straight line 
starting from the selected point is associated to each 
direction; (3) Nd equally spaced points are defined on 
each straight line. These points cover a certain maxi- 
mum distance, Lmax, measured from the selected cen- 
tral event. Each of these points -plus the fixed time tn 
is a possible reception event whose emission coordinates 
may be found by using the Newton-Raphson algorithm 
(see section 4) and, (4) from the emission coordinates 
(which allow us to find the satellite positions) we may 
determine the sign of x^ to know whether we are con- 
cerned either with a single (x^ < 0) or with a double 
(X^ > 0) positioning solution. 

For a given HEALPIx direction, the following cases 
have appeared: (i) x^ < from the chosen point on 
Earth {L = 0) to a certain distance L_ (where x^ van- 
ishes), and x^ > from L_ to L^ax, and (ii) x^ < 
from L = to L^ax, which means that either x^ does 



not vanish along this direction or Z/_ > L-max- We may 
then represent length L_ on a pixelised sphere by us- 
ing both the moUview projection (see Sect. 5) and an 
appropriate color bar; thus, the pixel color gives the dis- 
tance i_ for the corresponding direction. In this way, 
whatever the chosen event on Earth may be, we are 
displaying a surface which separates single from double 
positioning solutions for the fixed hypersurface = tn- 
Single solutions are located either on the mentioned sur- 
face or at the same side -with respect to the surface- as 
the Earth point chosen as a center. Double positioning 
solutions are all located at the other side of the surface, 
which is hereafter called the separating surface. 

In some papers (Coll, Ferrando fc Morales-Lladosa 
2010, 2012), the parts of the emission coordinate re- 
gion characterized by the conditions x^ < 0, x^ = 0; 
and x^ > 0, are denoted Cg, Ci, and Ct, respectively. 
Hence, our separating surface, the region containing the 
Earth point playing the role of a center, and the com- 
plementary region with double solutions are the inter- 
sections of the hypersurface = tn with Ci, Cg, and 
Ct, respectively. 

Figure 1 shows the separating surfaces corresponding 
to three 4-tuples of satellites belonging to the Galileo 
constellation. The value of t r has been arbitrarily cho- 
sen for each 4-tuple. Figure 2 displays the separating 
surfaces for the 4-tuple 1 considered in the top panel 
of Fig. 1. Six new values of Ir (six hypersurfaces) 
have been considered in this Figure (one in each panel) . 
They cover -together with the Ir value of the top panel 
of Fig. 1- a period of the Galileo constellation. In both 
Figures, the maximum distance has been chosen to be 
Lmax = 10^ Km, the white pixels contain the direc- 
tions pointing towards the four satellites at the cho- 
sen emission times, the garnet region corresponds to 
directions with L_ > Lmax, and the remaining pixels 
(Z/_ < Lmax) are colored according to the color bar 
appearing in each panel (numbers in the bar are values 
of L- given in Kilometers). It is easily observed that: 
(i) all the satellites seem to be included in a blue re- 
gion of influence, (ii) various satellites may be located 
in the same region, and (iii) inside the blue zones, the 
satellites may be located in the central part as well as 
in the zones close to garnet regions. Of course, the dis- 
tribution and positions of the blue regions around the 
satellites depend on the relative positions of the four 
satellites among them and with respect to the receiver 
(at emission times). In the directions of the garnet re- 
gions, our study has been repeated for a larger distance 
Lmax = 3 X 10^ Km, thus it has been verified that only 
for a very small number of pixels (located close to the 
non garnet region), a L- value satisfying the inequali- 
ties 10® < Z/_ < 3 X 10® has been found. Hence, in new 



6 



Figures corresponding to L^ax = 3 x 10^ Km, these 
pixels would not be anymore in the garnet region, but 
in the complementary one. Note that the new L^ax 
is close to the distance from Earth to moon and, con- 
sequently, if the moon is (is not) located in the garnet 
region of the chosen satellites, positioning on its surface 
would be single (double). In the case of double posi- 
tioning, additional measurement would be necessary to 
choose one of the two possible localizations. Positioning 
at these large distances -from the satellites- is theoret- 
ically possible, although technical problems would arise 
(weak signals, large positioning errors due to uncertain- 
ties in the satellite trajectories, and so on). 

Various sections, a;^ = constant, of some emission 
coordinate regions are represented in Figs. 1 and 2. 
Let us now consider the intersections of the celestial 
spheres of these Figures with planes containing meridi- 
ans; thus, 2D sections of the emission coordinate regions 
are found. It may be easily verified that, the structure 
of the resulting 2D sections is analogous to that of the 
sketch displayed in Coll, Ferrando & Morales-Lladosa 
(2012), but as expected, our realistic sections arc less 
symmetric. In both cases we may distinguish the 2D 
sections of the so-called central region with single po- 
sitioning (Cs U C|), and complementary 2D sections of 
Ct, which contain all the double locations. 

From Figs. 1 and 2, it follows that close enough to 
Earth, up to distance of the order of 10^ Km from the 
surface, positioning is single (the exact value of i_ de- 
pends on direction). However, for larger distances, the 
positioning may be either single or double and, con- 
sequently, if a Galileo (GPS) satellite is positioned by 
using four GPS (Galileo) satellites, intervals of single 
and double positioning are expected (see below). 

5.2 Co-region structure 

In order to study the co-region, four steps arc followed: 
(1) a point in the grid with coordinates (r^ ,t^,t^, tr) - 
which are the emission coordinates of the event selected 
in the last subsection (x^ < <;=^ single positioning)- 
is chosen to play the role of central point, (2) time 
r'* = Tr is fixed (section of the 4D co-region), (3) times 
T^, T^, and are varied along straight lines starting 
from the central point and following the 12288 direc- 
tions corresponding to a HEALPIx pixelisation (see 
Sect. 5), and (4) each line is uniformly covered by a 
set of Nd points, and quantity is calculated at each 
point to classify the associate reception event (single or 
double positioning). 

Points on the straight lines are characterized by the 
parameter A = [(r^ - rlf + (r^ - rlf + (r^ - rlfY'^. 
The values of this parameter are given in seconds. 



Along each straight line, it may be numerically verified 
that the emission-reception conditions (4) are only sat- 
isfied from A = to a cc^rtain value Xmax- Once qiian- 
tity \max has been numerically determined for each 
HEALPIx direction, the number of points and the 
separation between two neighboring ones may be appro- 
priately chosen to cover the segment limited by points 
A = and Xmax for any direction. 

We have verified that, along any HEALPIx direc- 
tion, quantity x^ is negative from A = to a certain A_ 
(where it vanishes), and positive from A_ to \max- Nev- 
ertheless, in the interval (A_,A„ia3;), where the condi- 
tion x^ > is satisfied, two cases may be distinguished: 
(i) there are two positioning solutions, and (ii) there are 
two future-like solutions (no positioning ones). 

In order to display the properties of the co-region, 
two panels are presented for each = tr hypersur- 
face. An appropriate HEALPIx pixelisation, the moll- 
wide projection, and color bars are used in both panels. 
In one of these panels the quantity represented in each 
pixel is A_ . Hence, we are representing a surface which 
surrounds the zone where the positioning solution is 
unique. This zone is hereafter referred as to the single 
valued co-region. In the other panel, the color bar is 
used to show the difference \max — A_ for the pixels cor- 
responding to directions with two positioning solutions, 
whereas the gray zone displays the pixels where no po- 
sitioning solutions exist for A > A_. In other words, the 
non gray part gives, for each pixel, the width of the zone 
where there are double positioning solutions, which is 
located outside the single valued co-region. This second 
external zone is hereafter called the double valued co- 
region. Fig. 3 shows the single (left panels) and double 
(right panels) valued co-regions for the same satellite 
4-tuples as in Fig. 1 and for arbitrary values of Ir. In 
Fig 4 (5) we display single (double) valued co-regions for 
the 4-tuple 1 and various times. The time tR of the top 
panel of Fig. 3 and the six times appearing in the panels 
of Figs 4 and 5 (the same considered in Fig. 2), cover a 
period of the Galileo satellites. The central point in the 
co-region and the equation of the hypersurface t"* = tr 
are obtained from the inertial coordinates of the chosen 
event on Earth by using the Newton-Rhaphson method. 
Since coordinates {x\,x1,x1) are the same in all cases 
(see above), only coordinate tR is given on top of each 
panel. In all cases, we find an internal single valued re- 
gion partially surrounded by a external double valued 
one. 

In the study of the emission region, the Newton- 
Raphson method must be applied for each direction of 
the HEALPIx pixelisation, whereas in the co-region, 
this method is only used at emission time. This makes 
the study of the region more time consuming. By this 
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reason, the chosen HEALPIx realizations of the region 
have less directions than those of the co-region (see 
Sees. 5, 5.1 and 5.2). Anyway, all the maps have a 
good enough angular resolution. 

5.3 Positioning GPS (Galileo) satellites with the 
Galileo (GPS) constellation 

It has been argued (Tarantola ct al. 2009; Coll, Fcr- 
rando & Morales 2006b, c; Coll, Ferrando & Morales- 
Lladosa 2010) that, in order to define an intrinsic co- 
ordinate system by using the proper satellites (with 
no reference to Earth), and also to measure the grav- 
itational field in the region where the satellites move 
(gravimctry), it is necessary the interchange of infor- 
mation among these satellites and with the receiver. 
Let us consider the most simple information exchange, 
in which one GPS (Galileo) satellite is considered as 
a receiver to be positioned with the emission coordi- 
nates broadcast by four satellites of the Galileo (GPS) 
constellation. On account of our study of the emission 
region presented in section 5.1 (see the last paragraph), 
it is evident that, in most cases, double positioning is 
expected to appear on the world line of the satellite 
playing the receiver role. 

We have first considered two cases in which a Galileo 
satellite (receiver) is positioned by using emission coor- 
dinates from four GPS satellites (emitters). In each 
case, the world line of the receiver is known. On ac- 
count of this fact, the following steps allow us to get 
the single and double receiver positions: (1) given the 
inertial coordinates of a point on the receiver world line, 
the Newton-Raphson method (see section 4) is used 
to get the emission coordinates, (2) from the result- 
ing emission coordinates, quantity is computed to 
determine the positioning character (single or double 
location) and, (3) if the positioning is double, we use 
Eqs. (5) to get both positions, one of them is always on 
the circular orbit of the receiver satellite, and the other 
one is outside this trajectory. Results are presented in 
Figs. 6 and 7. In order to build up these Figures, a 
special procedure has been designed which allows us to 
display the motion of both the receiver satellite along 
its circular orbit and the associated point (if it exists) 
on its trajectory. These Figures are 3D representations 
of the trajectories in the {x,y,z) = {x^,x'^,x^) space. 
Nc = 7200 equally spaced points are considered on the 
satellite circumference. Colors are used to follow the 
motion of the positioning solutions along their trajecto- 
ries. In Figs. 8 and 9, the same techniques arc used, and 
the same representations are displayed, but these Fig- 
ures correspond to a GPS receiver satellite positioned 
from the emission coordinates broadcast by four Galileo 
satellites. 



In Figs. 6 to 9, the motion along any trajectory is 
dextrogyre. Single positioning solutions are represented 
by red points; hence, these points are always on the cir- 
cular orbit of the receiver satellite. An initial point is 
arbitrarily chosen. Four subgroups of points have 
been selected. The double positioning solutions are dis- 
played by using the black, fuchsia, dark blue, and light 
blue, in the first, second, third and fourth subgroups, 
respectively. The first point of the first subgroup is the 
chosen initial point, and points and subgroups are or- 
dered in the sense of growing times (dextrogyre sense in 
Figures). According to these criteria, the initial point 
is marked by a red star on the circular trajectory for 
single positioning (as in Fig. 6 and 8), and by two black 
stars for double positioning (as in Fig. 7 and 9). The 
final point coincides with the initial one on the circum- 
ference, and the associated point (if it exists) will be ev- 
idently represented by a light blue star (as in Figs. 6, 7, 
and 9). Since the periods of the GPS and Galileo satel- 
lites are different, the character of the positioning (sin- 
gle or double valued) at the initial and final points may 
be different, it is due to the fact that, though these 
points coincide on the circumference, the locations of 
the four emitter satellites is different at the initial and 
final situation and, consequently, the positioning char- 
acter may be distinct. Moreover, as points on the cir- 
cumference -located inside a double valued region of 
any color- tend to possible points separating this re- 
gion of contiguous single valued ones, the associated 
point lying outside the circumference tends to infinity 
and, consequently, the two positioning solutions tend to 
a unique real one. Changes of color at the initial point 
appear as a result of differences between the positioning 
character at the initial and final situations (see above). 
In practice, the initial point never separates single and 
double valued zones and, consequently, if there is a fi- 
nal positioning point outside the circumference, it does 
not tends to infinite as the corresponding point on the 
circumference approaches the initial one. Taking into 
account all these criteria and comments. Figs. 6 to 9 
may be easily understood. The brief but illuminating 
description given in the Figure captions deserves atten- 
tion. 



6 Discussion and prospects 

This paper has been essentially devoted to the study 
of the bifurcation problem (double positioning) in rel- 
ativistic positioning systems. In order to develop this 
study, the following approach has been used: photons 
move in a 4D Minkowskian space-time, and satellites 
evolve in a 4D Schwarzschild space-time associated to 
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a spherically symmetric Earth. This procedure takes 
into account the effects of the Earth gravitational field 
on the clocks; for example, we have verified the well 
known fact that, for the GPS configuration, satellite 
clocks run more rapid than clocks at rest on Earth by 
about 38.4 microseconds per day. Our simulation of 
the GPS and Galileo constellations is accurate enough 
for many estimations and, in particular, to study bifur- 
cation. Double positioning situations arc identified on 
the emission coordinate region and co-region, as well as 
on the orbits of some satellites, which are considered 
as receivers to be localized with the help of four emit- 
ters. Methods to choose the true location in the case 
of bifurcation are proposed and discussed. Some open 
problems arc pointed out (see below). 

According to Coll, Ferrando & Morales-Lladosa 
(2010), the sign of is crucial to characterize single 
and double positioning. Since x is a vector orthog- 
onal to the hyperplane containing the four emission 
events, this hyperplane is time-like, null, and space-like 
for x^ > 0, x^ = 0, and x^ < 0, respectively. Single 
positioning corresponds to x^ < (space-like and null 
hyperplanes), whereas location bifurcations appears for 
X^ > (time- like hyperplanes). Once the emission co- 
ordinates are known, the internal satellite configuration 
and the sign of x^ may be easily obtained. Thus, the 
location bifurcations (double positioning, x^ > 0) may 
be found. In the double positioning cases, besides the 
emission coordinates, other measurements -either an- 
gles or times- have been proposed (see Sect. 3) to iden- 
tify the true location between the two possible ones. 
On account of these facts, we have found the zones 
in the emission coordinate region and co-region corre- 
sponding to single and double positioning. Results have 
been presented in Figs. 1 to 5. In any case, there is a 
central zone with x^ < O7 and a complementary one 
corresponding to double positioning. 

If a receiver (for example a satellite) is always located 
inside the central region of four emitter satellites, its 
positioning is single, and the emission coordinates are 
sufficient to find the receiver position; however, if the 
receiver enters and leaves the central zone either one 
or various times, there are phases of single and double 
positioning. These are the cases represented in Figs. 6 
to 9, in which, a satellite of a certain GNSS is positioned 
by using four satellites of other GNSS. A typical orbit 
of the GPS and Galileo satellite constellations enters 
and leaves various times the central region. This means 
that there are bifurcation locations and, consequently, 
apart from the emission coordinates, devices measuring 
either angles or times are needed to get the true receiver 
position. 

A clock aboard of the receiver satellite could mea- 
sure the observational time, to, when the emission co- 



ordinates are received. This measurement may be done 
whatever the positioning character may be. In the case 
of single positioning, there is only a coordinate time ti 
derived from the emission coordinates, whereas a pair 
of coordinate times (ti,f2) appears in any double po- 
sitioning. In the absence of errors in both the coordi- 
nate times and the observational time to, this last time 
would exactly coincide with ti for single positioning, 
and with one of the times (ti,t2) in the case of double 
positioning. The time coinciding with to would corre- 
spond to the true location. These coincidences would 
not be exact due to both the limited accuracy of the 
clock measuring to, and the positioning errors in ti or 
in the pair (ti,t2)- These errors may be associated, for 
example, to uncertainties in the world lines of the four 
emitters. The cases of single positioning might be used 
to calibrate the errors separating ti from to. For double 
positioning the differences \ti — to] and \t2 — to\ must 
be first estimated and, then, two cases may be distin- 
guished: (i) quantity \ti — ^2] is much larger than the 
typical value taken by |ti — fo| in single positioning, and 
(ii) the absolute value \ti —t2\ is of the order of the men- 
tioned typical value. In the first case, the smallest of 
the 1*1 —to\ and 1*2 —io| quantities clearly corresponds 
to the true positioning; however, in case (ii), the true 
position may not be found with the clock on board. In 
such a case, it could be studied if the criterion based 
on angle measurements (see Sect. 3) may lead to better 
results. 

We have estimated the maximum and minimum 
values of \ti — t2\ for the double positioning events 
appearing in Figs. 6 to 9. The maximum value is 
\ti—t2\max — 128.6 s. In general, large values of ^2! 
appear close to the transitions from double to single po- 
sitioning zones. There are no problems to choose the 
right position is these cases. The minimum value is 
1*1 — *2|mm — 1-1 X 10~^ s; it is also large as com- 
pared to the expected errors in the time to measured 
by a good atomic clock (aboard the receptor satellite). 
The errors due to uncertainties in the satellite orbits 
strongly depends on the Jacobian of the transforma- 
tion from emission to inertial coordinates (Puchades & 
Saez 2011). This Jacobian is to be calculated for the 
values of the emission coordinates corresponding to the 
double positioning under consideration. The detailed 
study of this type of errors as well as the analysis of 
other possible error sources are open problems requir- 
ing further research. 
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Fig. 1 All panels are HEALPIx moUwide representations of the distance L- in Kilometers. At this distance from the 
center, quantity vanishes for each direction in the hypersurface a:* = Ir. This surface separates single from double valued 
regions. The times tn and the 4-tuple of satellites used for positioning are given above each panel. 




Fig. 2 Same as in Fig. 1 for the 4-tuple 1 (top panel of that Figure) and for the tn times displayed above each panel. 
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4-tuple 1: tj;=13h 4-tuple 1; tR=13h 




Fig. 3 All panels are HEALPIx moUwide colored representations. Left (right) panels show the time distance A_ (Amax— A_) 
for each direction on the hypersurface — r* • The A_ values define a 2-surface where vanishes. Each point inside this 
surface leads to a single valued positioning in physical space-time. The gray zones of the right panels correspond to the 
directions for which there are no positioning solutions for A > A_. The colored part of these panels gives the width of the 
zone -surrounding the surface ~ 0- whose points lead to double valued positioning. Time tn and the 4-tuple of satellites 
used for positioning are given above each panel. 
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receiver tirQe=15h 
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Fig. 4 Same as in the left panels of Fig. 3 for the 4-tuple 



receiver time=17h 




receiver time = 21h 




0.041 0.24 
receiver time = 25h 




1 and the same tn times as in Fig. 2 



14 



receiver time=15h receiver time=171r 




Fig. 5 Same as in the right panels of Fig. 3 for the 4-tuple 1 and the same Ir times as in Fig. 2 




Fig. 6 Positioning a Galileo satelUte with four GPS emitters. The red star is a single valued initial point. It is the first 
point of a red arc of single valued solutions. Going through the circumference in the dextrogyre sense, we find a second 
continuous arc of double valued positions successively colored in fuchsia, dark blue and and light blue. This arc returns to 
the initial point. The line of the associated positions is overrun in the same sense. It tends to infinity in the transition from 
the fuchsia to the red arcs. The final positioning is double valued. One of the positions coincides with the initial one and 
the associated position is represented by the light blue star. 
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Fig. 7 Positioning a Galileo satellite with four GPS emitters. The two black stars are the initial points of a double 
valued solution. The following succession of arcs is observed on the circumference: black, red, fuchsia, red, light blue. The 
corresponding lines outside the circumference tend to infinity in the following transitions: from black to red, from fuchsia to 
red, and from light blue to red. As in Fig. 6, The light blue star is the final position. It is associated to the initial position 
on the circumference. 




Fig. 8 Positioning a GPS satellite with four Galileo emitters. The red star is the single valued initial position and also the 
final one. All the positions are single valued excepting a very small dark blue arc on the circumference. The corresponding 
curve outside the circumference tends to infinity in the two transitions from dark blue to red. 
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Fig. 9 Positioning a GPS satellite with four Galileo emitters. The two black stars correspond to a double valued initial 
position. On the circumference we find the following succession of arcs: black, red, fuchsia, dark blue, red, and light blue. 
The curves outside the circumference tend to infinity in the following transitions: from black to red, from fuchsia to red, 
from dark blue to red, and from light blue to red. The final point is represented by the light blue star associated to the 
initial position located on the circumference 



